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Abstract 



For Dirac neutrinos with magnetic moment, we compute the production 
rate for right-handed neutrinos in a hot and dense QED plasma containing an 
initial population of left-handed neutrinos thermally distributed. The most 
important mechanisms for vl depolarization, or production of right-handed 
neutrinos, are the — * vr chirality flip and the plasmon decay to &l + vr. 
The rates for these processes are computed in terms of a resummed photon 
propagator which consistently incorporates the background effects to leading 
order. Applying the results to the cases of supernovae core collapse and the 
primordial nucleosynthesis in the early universe, we obtain upper limits on 
the neutrino magnetic moment. 

PACS numbers: ll.10.Wx, 14.60.Lm, 95.30.Cq, 97.60.Bw 
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I. INTRODUCTION 



The properties of neutrinos have become the subject of an increasing research effort over 
the last years. Among these properties, the neutrino magnetic moment \i v has received 
attention in connection with various chirality flip processes that could have important con- 
sequences for the explanation of the solar neutrino problem |I]j2|], the dynamics of stellar 
collapse |3]|I|] and the evolution of the early universe ||. A non- vanishing neutrino mag- 
netic moment implies, for example, that left-handed neutrinos produced inside a supernova 
core during the collapse, could change their chirality becoming sterile with respect to the 
weak interaction. These sterile neutrinos would fly away from the star leaving essentially 
no energy to explain the observed luminosity of the supernova. The chirality flip could be 
caused by the interaction with an external magnetic field or by the scattering with charged 
fermions in the background, for instance vl e~ — > pr e~ and u^p — > urp. Invoking this last 
mechanism and using the average parameters inferred from the supernova 1987A, Barbieri 
and Mohapatra Q have derived a limit \l v < (0.2 — 0.8) x 10~ u fi B , where /xg is the Bohr 
magneton. 

It has also been pointed out that the constraints imposed by big-bang nucleosynthesis 
(BBN) do not allow the extra degree of freedom that a right-handed neutrino in equilibrium 
would introduce. In order to avoid that the chirality flip processes maintain a population 
of right-handed neutrinos in equilibrium during the evolution of the early universe, it is 
necessary that the average rate for these processes is less than the expansion rate of the 
universe at all times until the BBN epoch. Invoking this constraint, Elmfors et. al. || have 
derived a cosmological bound on the neutrino magnetic moment \i v < 6.2 x 10 _11 /xb. 

Dispersion processes in a plasma could exhibit infrared divergences due to the long- 
range electromagnetic interactions. To prevent such divergences, the authors in Ref. 0] 
introduced an ad hoc thermal mass into the vacuum photon propagator. However, it is well 
known that at high temperatures or densities, a consistent formalism developed by Braaten 
and Pisarski |7]|| and that renders gauge independent results, requires the use of effective 
propagators and vertices that resum the leading-temperature corrections. The method has 
been successfully applied to the study of the damping rates and energy losses of particles 
propagating through hot plasmas m this paper we use this framework to study the 

neutrino chirality flip processes in a dense and hot plasma. 

The most efficient process for conversion of left-handed to right-handed neutrinos hap- 
pens through scattering off electrons with the exchange of effective space-like photons. In the 
resummation method of Braaten and Pisarski, these photons are described by the spectral 
function of the photon propagator that develops a non-vanishing contribution for space-like 
momenta and whose physical origin is Landau damping. We compute the production rate 
of i/rs and the corresponding luminosity for such a process in a supernova. Our complete 
leading-order calculation is compared with the results obtained by means of an screening 
prescription used in a previous work f|. Our result can be used to place an upper bound 
on the neutrino magnetic moment which is in the range \l v < (0.1 — 0.4) x 10~ ll fi B [ 12 ]. 

This work is organized as follows: In section U we collect the ingredients that allow to 
compute the production rate of right-handed neutrinos from the imaginary part of the right- 
handed neutrino self-energy, in the real-time formulation of Thermal Field Theory (TFT), 
by means of a resummed photon propagator. In section Ell we restrict the analysis to the 
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production of right-handed neutrinos from the chirality-flip process. We obtain approximate 
expressions for the production rate that permit us to explore its analytical behavior in 
the small and large right-handed neutrino energy regions. In section IVl we study the 
production of right-handed neutrinos through the plasmon decay process which we show to 
be subdominant as compared to the chirality flip process. In section [V] we use the average 
parameters inferred from the supernova 1987A to find an upper bound to the neutrino 
magnetic moment. In section IVTL we also deduce an upper bound by imposing that the 
average production rate of right-handed neutrinos be at all times less than the Hubble rate 
up to the BBN epoch. This last result is shown to differ from that of Ref . j5| . We sumarize 
our results in section |V1J| and leave for the appendices some of the computations outlined 
throughout the rest of the work. 



II. FORMALISM 

Consider a QED plasma in thermal equilibrium at a temperature T such that T , jl e ^> 
m e , where m e and \x e are the electron mass and chemical potential, respectively. The pro- 
duction rate T of right-handed neutrinos with total energy E and momentum p can be 
conveniently expressed in terms of the z/r self-energy E as |T3" 



r(£)= ^^Tr[Pi?ImE] , (1) 

where L, R — | (1 ± 75) and rip is the Fermi-Dirac distribution for the right-handed neutrino. 
As shown below, in Eq. (0), the z/r thermal distribution cancels out from the final result. 
In what follows we consider E > corresponding to vr production. The annihilation of vr 
can be obtained from the case E < 0, however the initial vr population is negligible and the 
corresponding rate can be ignored. 

The resummation scheme is usually presented in the imaginary-time formalism. However, 
the expression for ImE can be directly computed following either the imaginary or the real- 
time formulations of TFT with identical results. In what follows we will work in the real-time 
formalism. As is well known, this formalism requires a doubling of the degrees of freedom 
and the propagators and self-energies adopt a 2 x 2 matrix structure. In particular, the 
imaginary part of the retarded self-energy is related to the 1-2 component of the self-energy 
matrix through 

ImE(P) = -^-E 12 (P), (2) 
lm F (E) 

where e(E) = 6(E) — 6(—E), with 9 the step function. As stated above the vr thermal 
distribution cancels out when Eq. (|2|) is substituted in Eq. ([]]). We take a = for the 
time path-parameter in the notation of Le Bellac || and Landsman and van Weert Jill • For 



simplicity we have selected the rest frame of the medium; however, the expressions can be 
rewritten in a covariant way replacing E by p ■ u, where is the velocity four- vector of the 
medium. 



3 




Fig. 1. Feynman diagram for the self energy E of the right-handed neutrino containing 
the effective photon propagator (resummed in the HTL approximation) denoted by a blob. 



The one loop contribution to S 12 , shown in Fig. 1, is given explicitly by 

Si 2 (P) = -inl J 0^ K *° aP S^(V+ H) LKpo** *D%(K). (3) 

We will use capital letters to denote four- vectors: P M = (E,p) for the incoming neutrino 
and = (ko, k) for the virtual photon, and p = \p\, k = \k\. For the neutrino-photon vertex 
we take the magnetic dipole interaction ^, v a a pK^ . The neutrino effective electromagnetic 
vertices are of no concern to us here since they are induced by the weak interaction of the 
particles in the background and thus conserve chirality ||15|| . 

The intermediate vl line can be taken as a bare fermion propagator because it gets 
dressed only through weak interactions with the particles in the medium. Hence, the S\2 
component for a massless fermion propagator is given by 

S 12 (Q) = 2m Q5 (Q 2 ) e(q ) n F (q ) , (4) 

where 

is the Fermi-Dirac distribution for the left-handed neutrino in the medium, with fi v being 
its chemical potential. 

In Eq. (^), the integration region where the momentum k flowing through the photon 
line is soft (i.e. of order eT) requires hard thermal loop (HTL) corrections to the photon 
propagator that contribute at leading order and must be resummed. The effective propagator 
is represented by the blob in Fig. 1, and is obtained by summing the geometric series of one- 
loop self-energy corrections proportional to e 2 T 2 . As usual, we split the photon propagator 
into longitudinal and transverse parts 

*D» V (K) = *D L {K) Pl v + *D T (K) , (6) 

we drop the term proportional to the gauge parameter since it does not contribute to E, as 
can be easily checked. In the previous equation all *D' 1U , *Dj- and are 2x2 matrices, 
while P£ v and P^ are the longitudinal and transverse projectors, respectively, 
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The effective photon propagator is obtained when the hard thermal loops in the photon 
self energy (first computed by Klimov [16] and Weldon |T7]]) are resummed. The complete 
matrix propagator is written as 



k D LtT (k) = U 











fc 2 -nt 



(8) 



where U is the photon thermal matrix |14J . The polarization functions IT^ and 11^ are given 
by 
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m 7 is the photon thermal mass, that in the limit T, jl e ^> m e is given by 

1 1 3#\ 
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(10) 



The Legendre function Qo(k /k) is defined in the complex k plane cut from —k to k; it is 
real in the time-like region but it acquires an imaginary part for space-like K 
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The solution of the equations K 2 — ReUi = and K 2 — ReUr = represent the propagation 
of longitudinal photons, or plasmons, and transverse photons, respectively. These solutions 
in the time-like region are well known, in general they are obtained numerically. However it 
is possible to obtain approximate analytical results for small and large values of k. In the 
small momentum limit k <C m 7 the dispersion relations reduce to 



<4 + \ k2 



(12) 



where the plasma frequency is defined as to p = J ~m 7 . In the large momentum limit k 3> m 
the behavior of the dispersion relations is approximated by 



UJrp = 171^ + k , 

ujl = k + 2fcexp 



~2 



m 



(13) 
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In a relativistic plasma the photon dispersion relation for the longitudinal mode K — 



RcUl = has also a solution in the space-like region fllSfl . Then the Cerenkov radiation of 
a plasmon is, in principle, kinematically allowed. The numerical solution of the longitudinal 
dispersion for u < k is shown in Fig. 2. We observe that the solution is close to the light 
cone. Using this fact and the first of Eqs. (Q), it is a simple task to derive the following 
approximate solution 

u>l = k — 2k -, , (14) 



1 + exp 



that, as shown in Fig. 2, agrees very well with the numerical solution. However, this mode 
develops a large imaginary part, which implies that the Landau damping mechanism acts 
to preclude its propagation. In this situation the correct method to include the complete 
contribution of both the space-like and time-like degrees of freedom requires the use of the 
spectral representation as given below in Eq. (|16|). The effect of the space-like mode is not 
very distinctive, except for the low energy spectrum of the Ur production via the spin flip 
reaction ui — > vr (see next section). 




Fig. 2. Dispersion relations for the longitudinal modes. The upper branch (dotted - 
dashed line) is the the usual time-like solution. The lower brach is the space-like mode. The 
numerical solution (solid line) is compared with the approximate analytical result (dashed 
line). 



From Eqs. @ and @ we can get the component *D 21 of the photon propagator. In 
order to obtain the complete result, we notice that for time-like momenta the II ^ and Ut 
functions are real, however for space-like momenta both functions acquire an imaginary part 
as seen from Eqs. and |TT 



U L.T 



The final result can be written as 
2m[l + f(k )]p LjT (k ,k), 



(15) 
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where f(k ) = (e fco//T — 1) 1 is the Bose-Einstein distribution and the functions Pl,t(^o ? k) 
are given by 



Pi (k ,k)= Zt (k) [S (k - LUi(k)) + 5(k + u^k))) + pi(k , k)6 (k 2 Q - k' 
for i — L,T. The residue Z^(k) for the longitudinal excitations is given by 

u L (k) 



Z L (k) 

and the cut function fl^ko , k) 
(5 L (k , k) = 



k 2 + 2m 2 



col(k) ' 
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(16) 
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where we have defined x = k /k. For the transverse mode, the residue and cut functions 
are given by 

u T {k) (u 2 (k) - k 2 ) 



Z T (k) 



(3 T (k , k) 



2m 2 u 2 T (k) - {u 2 {k) - k 2 ) 
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(19) 
(20) 



The functions Pl,t{<^ > k) are referred to as the photon spectral densities. It is an straight- 
forward exercise to verify that these coincide with the spectral densities computed in the 
imaginary-time formalism after the analytical continuation to real time is performed. The 
spectral densities Pl,t(u , k) contain the discontinuities of the photon propagator across the 
real-cu axis. Their support depends on the magnitude of the ratio between u and k. For 
\ui/k\ > 1, p LjT (u,k) have support on the points ±u LtT (k), i.e., the time-like quasiparticle 
poles. In the space-like region the support of Pl,t{^ k) lies on the whole interval —k < uo < k, 
with the contribution arising from the branch cut of Qq. Hence, the spectral density is the 
sum of pole and cut terms. 

To proceed with the calculation of Eqs. ([[]) and (|) we need the following quantities 

C T {E,K) = K a K p P TliV Tr [a a ^ (f + $) L f R] = k 2 (x 2 - if [{2E + u) 2 - k 2 



C L (E, K) = K a K p P Ltlu Tr \a a »(P +$)La^ PR 



k 2 (x 2 -l) (2E + u 



(21) 



The production rate of right-handed neutrinos in Eq. ([!]) can be computed after the substi- 
tution of the expressions for the fermion, Eq. (0), and photon, Eq. (|15D, propagators into 
Eq. (H). Using also the spectral representations, Eq. (|16|) , our result is 

° k e(E + h)t{h) (1 + f(k )) n F (E + k ) 



T(E) 



E 



(27TY 



[C L (E, K) p L (k , k) + C T (E, K) p T (k , k)} . (22) 

As we have already mentioned we shall consider E > corresponding to the production of 
v R . However, depending on the signs of E + oj and to the v R production rate can be divided 
into the chirality flip process Vl —>■ vr mediated by a virtual photon and the plasmon decay 
process 7 — > vl v r- We shall analyze these two processes separately in the following sections. 
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III. THE v L -> v R CHIRALITY FLIP 

The contribution to the vl ~~ > v r chirality flip process is obtained from the general 
expression in Eq. (p2| ) if we set E + uj > as corresponds to the case of an incident v L . 
The angular integration over the direction of k is readily performed. Using the condition 
| cosG| < 1, where is the angle between the momenta of the incoming neutrino and the 
virtual photon, we obtain two kinematical restrictions, namely, \u\ < k and (k — u) < 2E. 
The first condition implies that for the chirality flip ul —>■ vr process the kinematically 
allowed region is restricted to \u\ < k and consequently the only contribution from the 
photon spectral density in Eq. (|T6| ) arises from the cuts. The rate of production of right- 
handed neutrinos from the ul — » vr flip can be computed after the substitution of the 
expressions for the fermion, Eq. (g), and photon, Eq. (|1|), propagators into Eq. (|3]), the 
result is 

ir^ roo r-k 

T(E) = / kdk du;9(2E + uj-k)(l + f(uj))n F (E + u) 

16nE 2 Jo J-k 

[C L (E, K) /3 L (u, k) + C T (E, K) (3 t (uj, k)} , (23) 

where we set k = uj. Both, longitudinal and transverse photons, contribute to this rate. 
Notice that using the restriction (k — uj) < 2E, the integrand in the previous equation can 
be proved to be positive definite. The rate F can be written as the sum of two contributions 
T = T e + r a , that correspond to the production of Ur through the emission or absorption 
of a virtual photon. T e comes from the interval < uj < k, whereas T a corresponds to the 
interval —k < uj < 0, as can be checked by means of the identity 1 + f{uS) + f{—oS) = and 
the substitution oj — > —u in this second interval. 

The results for the longitudinal and transverse contributions to T{E) as a function of E 
are shown in Fig. 3, for a selection of values characteristic of a supernova core (see section 0). 
The solution for the chirality flip production rate requires in general numerical integration. 
However, the physical interpretation of the results is made clearer with approximate ana- 
lytical solutions that can be derived both in the small and large energy limits. 

Let us first consider the region of small neutrino energy (E T,ji v ). From Fig. 3, 
we observe that the longitudinal contribution shows a pronounced peak at energies below 5 
MeV. The result is somewhat unexpected because in the study of similar processes, the effect 
of infrared logarithmic singularities shows up in the transverse contributions. For example, 
in the calculation of the lifetime of a fermion quasiparticles in a QED plasma at high 
temperature the soft photon contribution leads to an infrared divergence in the transverse 
component alone due to the vanishing of the magnetic mass. The problem was solved with 
the use of a non-perturbative method based on a generalization of the Bloch-Nordsieck 
approximation In the present case, there are enough powers of k coming from the 

vertex factors in Eq. ([3]) to render the transverse contribution finite in the infrared. We will 
show this explicitly in appendix A. On the other hand the longitudinal photon is known to be 
screened by the Debye mass in the uj = 0, k — > limit, yet our results lead to a pronounced 
peak at small E. As we shall see, this effect arises because the longitudinal spectral density, 
Eq. (|T%|), has a logarithmic divergence at the light cone. A qualitative way to understand 
this behavior stems from the observation that there also exists a space-like branch for the 
longitudinal mode [IB], though with a large imaginary part that precludes propagation of the 
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mode. This branch lies very close to the light-cone. Could we ignore the large damping that 
the mode experiences, the situation would correspond to a process in which a left-handed 
neutrino changes its chirality becoming a right-handed neutrino through the emission or 
absorption of a Cerenkov photon. As we cannot ignore damping, the contribution shows up 
as that from a resonance that peaks for small right-handed neutrino energies, that is, for 
the region where the difference k — uj > is small. 




E (MeV) 

Fig. 3. The Longitudinal (solid line) and transverse contributions (dashed line) to the 
production rate T(E) for right-handed neutrinos via the spin flip transition —>■ Vr. 



Let us consider the limit E — ► 0, then the virtual photon is forced to be near the light cone 
(uj < k). The integration in Eq. (|23| ) is restricted to the kinematical region u < k < 2E + u, 
hence in the limit E — > we may set k = E + u everywhere in the integration. Keeping 
only the leading term in E, the cut contributions to the spectral densities can be replaced 
by the approximate expressions 

E 

(3t 



9 ' 



Pl = o , (24) 

2£m2 [In (£/T)f V ; 

where in the last expression we have simplified the approximation by replacing uj by its 
mean value T. We notice that the dependence on the longitudinal spectral function that 
originates the strong rise at small E comes from the fact that near the the light cone (3l 
behaves as 1/ [In (k — uj)] 2 . Using the previous results and the identity 

[1 + f(u)] [n F (E + u)] = n F (E) [f(uj) + n F (E + uj)] , (25) 

the resulting integrals can be evaluated analytically. Separating T(E) into its transverse and 
longitudinal parts, we obtain for the transverse contribution in the limit E —>■ 
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a 2 E 3 T 2 
Tt(E) m - n F (0) 
Tim* 



Li 2 (-e^ T ) - 



7T 2 



(26) 
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where Li n (z) = PolyLog[n, z] is the PolyLog function f2H |. For the longitudinal part the 
result reads 

T L (E) « ^ ^ (0) 2 [Li 5 (-e^) - C(5)l , (27) 

where is the Riemann zeta function. Both contributions vanish at E — 0, however the 
longitudinal part shows a steep rise for small E that originates the peak observed in Fig. 3. 
The maximum for the longitudinal spectrum can also be estimated and it is attained for a 
very small value of the energy 



E « T exp . (28) 




The previous results show that the Ul —>■ vr chirality flip reaction produces a huge number 
of small energy right-handed neutrinos. While this effect looks remarkable, the phenomena 
that it can give rise to are probably unobservable because at very low neutrino energies, 
experimental discrimination between the charged current interactions of low energy vr and 
vl would be fairly difficult. We also notice that in this limit, the reaction rate T(E) depends 
nonlinearly on the coupling constant e. 

We now turn our attention to the limit in which the neutrino energies are larger or of 
order T, jx v . As we shall see, a very good analytical approximation can be obtained in this 
case. The method extracts explicitly the leading logarithmic screening terms that arise from 
the use of the full resummed propagator. Following Braaten and Yuan |TTJ we introduce 



and intermediate cut-off q* such that m 7 <C q* T,jl v . In the region of hard momentum 
transfer k > q* the tree-level approximation is used for the virtual photon, q* acting as an 
infrared regulator. In the soft region k < q*, the effective resummed propagator is used. 
Adding the hard and soft contributions, the dependence on the arbitrary scale q* cancels. 

Let us first outline the calculation of the soft-momentum transfer contribution to T{E). 
In this region, hard thermal loop corrections to the photon propagator are not suppressed by 
powers of e, hence the resummed photon propagator must be used. Restricting our attention 
to the momentum region k < q* we can set 1 + f(ui) ~ T/u and for the functions in Eq. (|2~T|) 
we make the approximations 

AK A E 2 

C T (E, K) « -C L (E, K) « — — (29) 
inserting the results into Eq. (|23"D, this last reduces to 

v, ol m = H^SL £ k * dk £ £ (i _ 2 k) _ k)] . (so) 

The integral over u can be evaluated by using the sum rules derived from the analytical 
properties of the effective propagators |21| 
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The logarithmic dependence on q* can be extracted analytically, the result is 
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(32) 



(33) 



where we used the fact that m 7 <C q* and have extended the limit of integration to oo since 
we are just interested in the leading contribution. The remaining integral inside the brackets 
is a dimensionless quantity and can be computed numerically with the result —0.14. 

We now turn to the hard- momentum transfer contribution to T(E). In this region the 
tree-level photon propagator must be used. Following Le Bellac || we notice that tree-level 
results can be recovered from the full calculation by simply neglecting m 2 in the denomina- 
tors of Eqs. (^) and (|2T)D. Hence, the approximated spectral densities yield 



(3 L « 2{3 T 



m^ui/k^ 
k 2 — uj 2 



(34) 



It is convenient to further decompose the hard region (k > q*) into a (I) low {u < q*) 
and (II) high (uj > q*) frequency regions. In the low frequency region we can still use 
the approximations 1 + f{uj) ~ T/u and Eq. (|2"5|). The remaining integrations are readily 
performed. Recalling that q* E, the result is 



1 hard(E) ~ ^ n F (E) . 

37T 



(35) 



In the high frequency region it is not obvious that we can neglect u and k as compared to E, 
in particular the approximation in Eq. (|29|) may not be accurate. However for E ~ T,fl u the 
integrals are dominated by small uj and k, hence neglecting u>, k as compared with E yields 
the the leading order contribution in T/E. The approximation can be improved for large 
E, by computing the next order term, the second order corrections are given in appendix 
B. Thus to leading order we can set u>, k — > as compared with E wherever possible, except 
for the Bose and Fermi distributions that provide the cutoff for the integrals and must be 
taken exact. Using the identity in Eq. (E5f) we obtain 



r ii 
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(36) 



The complete rate for the production of right handed neutrinos to leading order in e is the 
sum of Eqs. (BO), (0) and fl35D. The result is 
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T{E) = 2?r r n F {E) 



C + ln 



1 / cosh 2 sinh ( 

+ 2 ln { cosh (J) 




(37) 



where C = 1.88. Note that the dependence on the arbitrary intermediate scale q* cancels 
out between the soft contribution in Eq. ( |30D and the hard contribution in Eq. (|35|). In 
Fig. 4, we display the comparison of this analytical approximation with the result of the 



numerical integration of the exact result in Eq. (23). We notice that for energies in the 
region E 3> T,£i u the analytical result gives an excellent approximation to the exact result. 
As mentioned before the result can be improved by computing the next order correction in 
T/E. Fig. 4 also shows the result of the second order contribution calculated in appendix 
B. 




Fig. 4- The production rate T(E) for right-handed neutrinos via the spin flip transition 
v l v r- The exact (numerical) result in the solid line is compared with the approximate 
result (dashed line) to leading order in T/E and the approximate result (dotted- dashed line) 
to second order in T/E. 



The energy carried by the produced right-handed neutrinos or vr emissivity q can be 
obtained from T(E) according to the relation 

q= iw mE) - < 38) 

The Ur energy spectrum (oc E 3 T(E)) as a function of E is shown in Fig. 5, where the exact 
numerical result, shown as solid line, is compared with the approximate analytical results. 
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We observe an excellent agreement between the exact and the second order approximate 
solution. 

Finally, we quote that in the degenerate limit T — > 0, T(E) reduces to a very simple 
expression 



r(E) 



Air 



(fiv - E) 9(fi v - E) 



(39) 



which is similar to the one obtained in Ref. [^2| for the damping rate of a fermion in a QED 
dense relativistic plasma at zero temperature. 
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Fig. 5 The energy spectrum ~ E 3 T(E) for right-handed neutrinos produced by the spin 
flip transition ul — > vr. The exact (numerical ) result in the solid line is compared with the 
approximate result (dashed line) to leading order in T/E and the approximate result (dotted 
line) to second order in T/E. 



IV. PLASMON TO v L v R DECAY 

We now consider the contribution to the vr production rate arising from the plasmon 
decay process 7 — > vlVr- The rate is obtained from the general expression in Eq. (|22| ) by 
setting E + k Q < as corresponds to a negative energy ul or an outgoing v^. The angular 
integration over the direction of k leads now to the following kinematically restrictions 

uj = —k > 

uj- <E <uj + . (40) 

where — (u ± k)/2. These restrictions imply that for the plasmon decay process, the 
kinematically allowed region corresponds to the time-like region \u\ > k and thus the con- 
tribution from the photon spectral density in Eq. ([16]) arises solely from the poles. Using 
the identities 
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1 + f(-Cu) = -f{ U ) 

n F {-E>) = 1 - e{Et+ ^ )/T - i = 1 - n F {E>) , (41) 

where rip refers to the thermal distribution of we can obtain the vr production rate 
from the plasmon decay process. Substituting the expressions for the fermion and photon 
propagators, Eqs. (H) and (|15|), into Eq. ([|) we obtain 



r ( E ) = 77rk E rkdk9(u + -E)9(E-cu-)f(u; i (k))[l-h F (u; i (k)-E)} 

dfa-UiikfiZiik), (42) 

where the sum is carried over the longitudinal and transverse modes. The vr emissivity can 
now be computed using Eq. (p8|) . In general, the E and k integrations have to be computed 
numerically, however, in many applications we can consider that the r?i neutrinos are absent 
from the medium. Then, the Pauli blocking factor can be neglected, and the integration 
over the Vr energy can be readily performed with the result 



2 



967T 3 Jo 



OO 



k dk 



1 / N 2 



-u L (k) (u 2 L (k) - k 2 ) f(uj L (k))Z L (k) 
+u T (k) (u*(k) - k 2 ) 2 f(u T {k))Z T (k)\ . (43) 



The remaining integral still requires numerical computation. However, analytic expressions 
can be obtained in the limit cases m 7 <C T and m 7 ^> T. In order to discuss these limits, 
it is convenient to decompose the emissivity q into its longitudinal and transverse parts qL 
and qr- 



The limit m 7 T requires ultrarelativistic temperatures and, according to Eq. (1C), an 
electron density bound by the condition efi e \^2nT. In this limit, the integral for the 
transverse emissivity in Eq. ( |4"3]) is dominated by k on the order of T. Thus, we can use 
the large k limit of the dispersion relation in Eq. (|T3| ) and the resulting integral can be 
evaluated analytically with the result 

" T = w 1 " 1 ^ 3 - (44) 

where ((z) is the Riemann zeta function. 

For the longitudinal emissivity in the same limit m 7 T, the contributions for large 
k can be neglected. The reason is that the substitution of the asymptotic expansion, Eq. 
(fnj), in the factor (^(/c) 2 — k 2 ) 2 produces a Gaussian cutoff in the integrals. Thus, the 
integral is dominated by k on the order of m 7 or smaller. The Bose distribution can be 
approximated by /(^x) ~ T/uj^. The integral then involves only the scale m 7 and can be 
evaluated numerically 

qL = o - 28 S m ° T - (45) 

We notice that in this limit the longitudinal emissivity is negligible as compared to the 
transverse emissivity, since it is suppressed by a factor (m 7 /T) 2 . 
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The limit T <C m 7 is obtained for an ultrarelativistic electron density and low tempera- 
tures according to the relation e/2 e ^> \f2-KT. In this limit the integrals in both the longitu- 
dinal and transverse cases are dominated by momenta small compared to m 7 . We can set 
k —>■ wherever possible except that, according to the asymptotic expansion in Eq. (|12|). 
the Bose distribution is approximated by / — > exp[— {u p + \k 2 / uj p ) /T] with A = 3/5, 3/10 for 
the transverse and longitudinal modes respectively. The resulting integrals can be evaluated 
analytically with the results 

Q T = ^ (1^1) 1/2 w ll/2 T 3/2 -u, p /T 
qT 768vT 3 V 27 7 P 



uip/T 



V. THE v R EMISSION IN A SUPERNOVA 

As a first application of our results we consider the emission of right-handed neutrinos 
immediately after a supernova core collapse. The large mean free path of the right handed 
neutrinos compared to the core radius implies that the vrs would freely fly away from the 
supernova. Therefore, the core luminosity for vr emission can be simply computed as 

Qvn = Vq, (47) 

where V is the plasma volume and q is the vr emissivity computed from Eq. (|38|). To make 
a numerical estimate, we shall adopt a simplified picture of the inner core, corresponding to 
the the average parameters of SN1987A p3|j24| . Consequently, we take a constant density 
p « 8 x 10 14 g/cm 3 , a volume V ~ 8 x 10 18 cm 3 , an electron to baryon ratio Y e ~ Y p ~ 0.3, 
and temperatures in the range T = 30 ~ 60 MeV. This corresponds to a degenerate electron 
gas with a chemical potential jl e ranging from 307 to 280 MeV. For the left-handed neutrino 
we take pL v ~ 160 MeV. Using this values in Eqs. (|2"3"| ) and P7|), we obtain by numerical 
integration 

Q VR = (^j (0.7 - 4.3) x 10 76 ergs/sec , (48) 

for T ranging from 30 to 60 MeV. The main contribution to this result arises from the Ul — > 
vr flip process, whereas the contribution from the 7 — > v^vr decay is smaller by two orders 
of magnitude. This result is in agreement with the observation first made by Fukugita and 
Yazaki [551] who noticed that for cosmological and astrophysical scenarios the plasmon decay 
process is subdominant as compared to the chirality flip process. Moreover, an estimate 
of the vr luminosity derived from the approximate solutions discussed in section [III] are 
surprisingly accurate, they differ from the result in Eq. fl48|) by less than 2%. 

Assuming that the emission of vrS lasts approximately for 1 sec, the luminosity bound 
is Q VR < 10 53 ergs/sec. which places the upper limit on the neutrino magnetic moment 

\l v < (0.1 -0.4) x lQ- U fx B . (49) 
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This upper bound slightly improves the result previously obtained by Barbieri and Mohap- 
atra ||. As mentioned before, these authors consider the helicity flip scattering u L e — > v R t 
to order e 4 introducing the Debye mass in the photon propagator as an infrared regulator. 

A word of caution should be mentioned in relation to the result in Eq. (|49|). It has been 
pointed out by Voloshin |26| that the z/r's produced by the magnetic moment interaction 
could undergo resonant conversion back into u^s through spin rotation in the magnetic field 
of the supernova core, with the subsequent trapping of the u^s by the external layers. If 
this is the case, then the bound in Eq. (f49[) becomes meaningless. However, the core density 
is rather high and the matter effect might dominate over the \i v B term, suppressing the flip 
back of v R to v L [p3f| . 

Recently, another mechanism for the neutrino chirality flip has been proposed, which 
occurs via the Cerenkov emission or absorption of plasmons in the supernova core p? 



Since the photon dispersion relation in a relativistic plasma shows a space-like branch for 
the longitudinal mode, the Cerenkov radiation of the plasmon is, in principle, kinematically 
allowed |Tj|. However, this mode develops a large imaginary part, which implies that the 



Landau damping mechanism acts to preclude its propagation as we have discussed. Conse- 
quently, we think that no better than the quoted limit in Eq. (|||) can be derived by this 
type of neutrino chirality flip processes in a supernova core. 



VI. EARLY UNIVERSE 

As a second application of our results, we now consider the production of right-handed 
neutrinos during the evolution of the early universe. If the rate of production of these 
neutrinos is able to maintain them in thermal contact with the rest of the particles in the 
plasma, then they will contribute to the effective number of degrees of freedom until their 
final decoupling. This could in principle affect primordial nucleosynthesis. In order to 
prevent right-handed neutrinos from being in thermal equilibrium, we need to require that 
their average production rate be less than the Hubble rate at all times until the neutrino 
freeze-out epoch. During the radiation dominated era, the Hubble rate was 




H = [ — g* ) , (50) 

mpiank 

with ~ 10.75 the effective number of degrees of freedom at the nucleosynthesis epoch. 
On the other hand, the average right-handed neutrino production rate can be obtained from 
Eq. (|23| ) (neglecting the contribution from the plasmon decay process) by averaging with an 
equilibrium distribution 

n u (E) = ) (51) 

exp{E/T) + 1 

appropriate for the early universe where the chemical potentials should be negligible. There- 
fore, the condition to avoid populating the right-handed neutrino component becomes 




<r> = 5.78 x 1(T 4 T 3 //* < — ( ] , (52) 
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where for nucleosynthesis, 100 MeV > T > 1 MeV. The most stringent bound on [i v is 
obtained for the highest possible temperature and thus we take T = 100 MeV which yields 

//„ < 2.9 x 10~ 1c Vb . (53) 

This result has to be compared to that from Ref. || where the use of a full one-loop ap- 
proximation to the photon polarization functions is used, instead of Eqs. (^). This choice 
leads to an upper bound one order of magnitude smaller than the above. One should notice 
however that when a plasma is such that the largest energy scale available is set by the 
temperature (or density), as in the present scenario, it is necessary to select the leading 
temperature (density) contributions out of perturbative calculations and sum these up into 
effective vertices and propagators in order to extract from them meaningful quantities. In 
this manner, one can ensure that the leading perturbative corrections are effectively taken 
into account [/J. Therefore, we conclude that the bound on the neutrino magnetic moment 
set by nucleosynthesis constraints is not nearly as stringent as that set by the analysis of a 
supernova core collapse. 



VII. CONCLUSIONS 



In conclusion, we have shown that the scattering processes mediated by effective plasma 
photons allows for the efficient conversion of vl into vr. In this work, plasma effects are 
consistently taken into account by means of the resummation method of Braaten and Pis- 
arski within the real-time formulation of TFT. For soft values of the energy, the production 
rate of vr differs significantly from that obtained by a constant Debye mass screening pre- 
scription. However, corrections to the integrated luminosity are small and for this reason, 
our upper bound on the neutrino magnetic moment does not differ significantly from the 
one obtained by Barbieri and Mohapatra H. Knowledge of an accurate expression for the 
Vr production rate, as given in Eq. (p3|), could be of importance in a detailed analysis of 
supernova processes. We also obtain another constraint on \i v by considering the possible 
effect of vr production in the early universe. The upper bound imposed by the analysis of 
SN1987A is two orders of magnitude smaller than the one obtained from the nucleosynthesis 
constraint. 



VIII. APPENDIX A 

In this appendix, we want to explicitly show that the transverse photon contribution, 
Tt{E), to the right-handed neutrino production rate, T(E), is free from infrared divergences. 
To this end, we refer back to Eqs. (p0|) and (p3|) that give the explicit expression for T T (E). 
Let us write x = u/k thus, when x — > 0, we can write 

l + /(*s)->^, (54) 

where f(z) is the photon statistical distribution. Notice that in this limit, /3x(k, x)[l + f(kx)} 
remains finite, unless we also take k — > 0. Let us write 
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i3 T (k,x) 



1 

k 2 (l-x 2 



-/3 T (k,x) 



with the definition 



k 2 mlx(\ - x 2 ) 2 /2 



k 2 {l - x 2 ) + m 2 (x 2 + f (1 - x 2 ) In (£5))" 
Notice that for x <C 1, 

(3 T (k,x) m 2 k 2 

* ~T k 4 + (nm 2 x/2) 2 ' 



+ 



2 (1-z 2 ) 

irm^x 2 ; 



Therefore, in the limit k — ► 0, 



(3r(k,x) 



x 



S(x). 



r soft 

1 T 1 



16tt£ 2 

16tt£ 2 



rip(E) I dkkO{2E -k)[{2E) 2 -k 2 
n F {E) 



o 

m 2 (2£ 2 -m 2 /4), 2£ > m 
4£; 4 , < m 



7 



7 



(55) 



(56) 



(57) 



(58) 



In this manner, we can write the contribution to Tt{E) from the soft (k — > 0, — > 0) region 

as 



(59) 



where in the integration we have set as the upper limit the soft scale m T This choice is 
somewhat arbitrary but it does not matter as long as for the hard contribution to Ft(E) we 
integrate k from m 7 . Eq. (|59D is thus explicitly infrared finite. 



IX. APPENDIX B 

The second order correction in (T/E) to T(E) for the quirality flip ul — > vr reaction 
can also be calculated analytically As discussed in section |T|, the leading order term is 
obtained by setting lu, k — > as compared with E wherever possible, except for the Bose 
and Fermi distributions. In the soft-momentum region, the condition uo,k E is valid 
everywhere, and consequently there are no T/E corrections to the result in Eq. (|55|). 

The hard-momentum region (k > q*) was decomposed into a (I) low (u < q*) and (II) 
high frequency (u > q*) regions. In the low frequency region we have to <C E and only k/E 
corrections have to be computed. It is a simply exercise to show that the only corrections 
are of order (T/E) 2 , hence the result in Eq. (|35|) is not modified in the next order. Finally, 
in the high frequency region we compute corrections of order k/E and lu/E to the result in 
Eq. (|36"|). In particular, instead of Eq. (p9|) , we approximate the functions in Eq. (|2l| ) by 

C T {E, K) w -C L {E, K)^^- (E 2 + Ecu) . (60) 
The correction to the result in Eq. ( |36D is given in terms of PolyLog functions as 
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T(E) = — - — — np{E) 



7 



2n 
E 
T 



In 



T 
8E 

1 + e^' T " 

1 - p-E/T 



) (-7T 2 + 7Li 2 (e~ E ' T ) + 7Li 2 (-e^ T ) - 8Lz 2 



7 (E\ 



2 n 



2 VT 



(61) 



This corrections improve the approximate analytical solution given in Eq. fl3"7|) when com- 
pared with the exact result in Eq. (|23|) for energies E > T, fj, u , as shown in Figs. 4 and 
5. 
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